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Collisional  Relaxation  of  Non-Equilibrium  Vibrational  Energy  Distributions 
in  a Thermal  Unimolecular  System.9  Surface  Collisions. 

by  D.  F.  Kelley,  B.  D.  Barton,  L.  Zalotaib  and  B.  S.  Rabinovltch 

Department  of  Chemistry,  University  of  Washington,  Seattle,  WA  98195 


There  have  been  no  measurements  to  date  of  the  dynamics  of  collisional 
relaxation  of  vibrational  distributions  characteristic  of  reacting  polyatomic 
molecules  at  high  energies.  Formal  theoretical  treatments  have  been  given. ^ 
The  unknown  desideratum  has  been  the  appropriate  matrix  of  the  collisional 
transition  probabilities,  P.  The  same  matter  appears  In  the  form  of  the  mean 

first  passage  time,  x , of  molecules  with  respect  to  an  absorbing  barrier 

3 4 — 

at  energy,  EQ.  Kim  and  Widom  have  given  expressions  for  x in  terms 

of  P . For  a thermal  unimolecular  system,  the  concept  of  an  absorbing  barrier 

corresponds  to  the  limiting  low  pressure  behavior.  For  the  trivial  case  of 

a constant  temperature,  steady-state  system,  x = l/kQ , where  kQ  is  the  low 

pressure  unimolecular  rate  constant.  But  x has  not  been  measured  for  the 

dynamical  case  of  transients  wherein  some  initial  input  distribution,  not 

necessarily  Boltzmann,  relaxes  to  a final  distribution  characteristic  of  the 

bath  temperature,  whether  by  homogeneous  or  heterogeneous  collision  processes. 

5 

Based  on  a quasi -accommodation  model  for  energy  transfer,  as  well  as  experi- 
mental results  with  bath  molecules  of  Increasing  complexity,®  one  anticipates 
that  the  wall  may  be  at  least  as  efficient  as  any  bath  molecule. 

A very  simple  technique  for  studying  the  dynamics  of  energy  transfer  at 
a surface  has  been  developed.7  The  apparatus  consisted  simply  of  a large, 
degassed  5 £.  fused  quartz  bulb,  modified  with  several  attached,  heated 
"reactor"  zones  of  varying  dimensions  (cylindrically  symmetric)  distributed 
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over  the  surface.  The  bulb  was  connected  to  a gas  handling  system  and  could 
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also  be  heated  to  alter  the  input  distribution.  Entrance  area  to  a reactor  was 

restricted  to  < 1/100  of  the  bulb  area.  Dimensions  of  the  reactors  used  so  far 

provided  mean  wall  collision  numbers,  m,  of  2,  4.5,  10.5  and  22.3 , per  encounter. 

For  a given  m , the  distribution  of  wall  collisions  was  found  by  Monte  Carlo 
7 8 

calculation.  * Molecules  that  returned  to  the  bulb  were  cooled  by  collisions 

before  another  encounter  with  the  reactor.  Any  value  of  in  from  1 to  10^  Is 

readily  attained;  the  technique  is  termed  the  Variable  Encounter  Method. 

We  have  applied  the  method  to  the  structural  Isomerization  of  1,1-cyclo- 

(E0  a 63  kcal  mole-1). 

propane-d2  activated  by  collisions  at  a seasoned  surfaceA  This  system  is  notably 
free  from  complications.  Brief  pre-exposure  of  the  oxidized  silica  wall  to  sub- 
strate provided  adequate  seasoning  and  reproducible  results.  To  perform  a run, 
one  reactor  was  heated  to  a desired  temperature,  the  bulb  (at  375  K)  was  pumped 

C 

to  a stable  residual  pressure  of  < 10  torr,  and  substrate  was  admitted  at  a 
-4 

pressure  of  1-2x10  torr  and  left  for  a suitable  time.  A flow  mode  was  also 
possible.  Analysis  was  by  gas  chromatography;  percent  reaction  varied  from 
1 - 50%. 

Figure  1 shows  experimental  temperature  plots  of  the  average  reaction  proba- 
bility per  collision,  Fc  , for  the  four  reactors.  The  data  were  compared  to 
stochastic  calculations  using  two  forms  of  P , gaussian  and  exponential.  Several 
models  were  computed  In  which  the  average  down jump  step  size,  <AE>  , varied  as 

Ex,  x * 0,  1/2,  1.  The  gaussian  form  (with  x = 0,  1/2) tentatively  give  the  best 
(cf.  refs.  6 and  9). 

1 f1tsA  Sorie  numerical  values  of  <AE>  for  the  case  x 3 0,  are  given  In  the  figure 

(use  of  x = 1/2,1  leads,  of  course,  to  larger  values  of  <AE>  at  the  level  of 

, _ 34 

energy  EQ) . Once  £ has  been  selected,  t follows  Immediately.  * 

Figure  2a  Illustrates  the  variation  of  Pc  with  m at  constant  temperature. 

Y 

Figure  2b  shows  the  sequential  variation  of  Pc  vs  collision  number,  calculated 
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from  the  flat  gaussian  model. 

We  find:  a)  For  the  temperatures  covered,  accommodation  to  a steady 
state  with  the  wall  takes  place  in  a small  number  of  collisions  % 8-15  (Fig.  2b), 
the  number  increasing  with  temperature;  this  accounts  for  the  findings  made  in 
low  pressure  pyrolysis  studies10  that  m values  £ 80-100  provide  temperature 
accommodation;  at  lower  temperatures,  especially,  wall  efficiency  approaches 
strong  collider  behavior  and  t is  quasi -independent  of  the  Input  distribution 
(bulb  temperature),  b)  Wall  collisions  are  more  efficient  than  gas  phase 
collisions;  the  surface  seasoned  with  a film  of  "polycyclopropane"  Is  charac- 
terized by  down  step  sizes  <AE>  of  4,000  and  3,400  cm"1,  (x  = 0),  at  775  K and 
11 

975  K,  respectively;  these  compare  with  the  lesser  values,  ^ 3,500  and  1,900 
-1  12 

cm  , measured  previously  for  binary  gas  phase  cyclopropane-d^  collisions; 

c)  The  present  data  thus  support  the  previous  finding  that  energy  transfer 

collision  efficiency  declines  at  higher  temperatures  for  this  polyatomic  uni- 

1 3 

molecular  system.  Dove  and  co-workers  have  estimated  an  "Incubation"  period 
of  'v  3,000  collisions  for  the  shock-heated  trl atomic  NgO  system  above  2000K. 

Features  such  as  surface  structure,  substrate  structure  and  level  of  EQ 
are  subjects  for  future  investigation  by  this  technique. 
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Figure  Captions 


1.  Plots  of  experimental  values  of  ?c  vs  T for  each  of  the  four  reactors. 
Solid  curves  are  calculated  fits  for  a gaussian  model  (x  = 0)  of  the 
energy  transfer  probabilities  for  each  value  of  m,  and  include  the 
limiting  case,  m = °°  . The  best-fit  values  of  <AE>  (cm~^)  are  given 
at  several  temperatures. 

2.  A)  Illustrative  plots  of  experimental  values  of  Pc  vs  m at  823K,  973K 
and  1123K.  Solid  curves  are  calculated  for  respective  values  of  <AE> 
of  3800,  3400  and  3050  cm"^.  The  ordinate  scales  are  arbitrary  and  are 
normalized  to  'P  = 1 for  m = °°  , at  each  temperature. 

B)  Illustrative  histograms  of  the  calculated  sequential  probability 
Pc(n)  v^  n , the  number  of  consecutive  collisions,  for  the  temperatures 
and  values  of  <AE>  used  in  A). 
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